arXiv:1504.04570v2 [cond-mat.supr-con] 14 Jan 2016 
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In a fully-gapped superconductor the electronic Raman response has a pair-breaking peak at 
twice the superconducting gap A, if the Bogoliubov excitations are uncorrelated. Motivated by 
the iron based superconductors, we study how this peak is modihed if the superconducting phase 
hosts a nematic-structural quantum critical point. We show that, upon approaching this point by 
tuning, e.g., doping, the growth of nematic correlations between the quasiparticles transforms the 
pair-breaking peak into a nematic resonance. The mode energy is below 2A, and stays hnite at the 
quantum critical point, where its spectral weight is sharply enhanced. The latter is consistent with 
recent experiments on electron-doped iron based superconductors and provides direct evidence of 
nematic correlations in their superconducting phases. 


Introduction .— Nematicity in correlated metals im¬ 
plies the propensity of the interacting electron fluid to 
break rotational symmetry without necessarily breaking 
translation symmetry [T]. It is at the heart of the phase 
competition in iron based superconductors (Fe SC) [2], 
leading to a lattice softening EHa, transport and mag¬ 
netic anisotropies iHin], signatures in optical [TTl IT2] . 
single particle mHis] and Raman m spectroscopies in 
the normal state. In the latter, the critical nematic fluc¬ 
tuations in the tetragonal phase appear in the Big Ra¬ 
man response as a low energy quasi-elastic peak. 

While nematicity in the normal state of Fe SC has 
been studied extensively, it is unclear whether nematic 
fluctuations couple to low energy charge carriers in the 
superconducting state. In fact, even the existence and 
detailed nature of the nematic fluctuations below Tc is 
not established either experimentally or theoretically. 

In this Letter we demonstrate that, near a nematic 
quantum critical point (QCP) which is well inside a fully- 
gapped superconducting phase [IZIIIH], the growth of ne¬ 
matic correlation transforms the standard pair-breaking 
peak in the Raman Big channel into a new nematic 
resonance mode at finite frequency. In a clean system 
the resonance is a sharp structure below 2A of the elec¬ 
tron pockets of Fe SC, and its spectral weight increases 
strongly near the QCP. The latter, which is a key feature 
of the theory, and which cannot be understood within 
a pair-breaking peak scenario, is consistent with exist¬ 
ing Raman scattering measurements in Co doped Ba-I22 
[191 [20] Na-111 [21] systems. This gives very strong 

evidence for this new mode, and, in turn, reveals that 
there is a substantial coupling between Bogoliubov quasi¬ 
particles and nematic fluctuations. 

When compared to the normal state behavior, the res¬ 
onance in the superconducting phase can be understood 
as the consequence of the opening of the gap, which shifts 
spectral weight to higher frequencies, thereby transform¬ 
ing the quasi-elastic peak into a resonance. It is caused 


exclusively by electronic contributions, while renormal¬ 
izations due to the lattice are dynamically screened and 
are thus suppressed. This leads to the interesting result 
that the Raman nematic resonance energy remains finite 
even at a nematic QCP. 

The nematic resonance is analogous to the spin- 
resonance, observed in neutron scattering experiments 
in Fe SC and cuprates superconductors [22ll27] . In fact, 
there is a close formal analogy of our theory to the one 
for the spin-resonance mode in superconductors with sign 
changing gap near a magnetic instability [28l[29]. In both 
cases a singular behavior of the imaginary part of a bare 
response function, caused by the superconducting coher¬ 
ence factors, is amplified by the vicinity to a QCP. Im¬ 
portant differences are that the singularity in our case 
is not related to a sign change of the gap as small mo¬ 
mentum excitations are probed. For the same reason the 
Raman resonance is more sensitive to nodes of the gap. 

Nematic resonance .— In the following we develop a 
theory for the nematic resonance, and we argue that it 
is intimately related to the quasi-elastic peak with di¬ 
verging spectral weight experimentally observed in the 
normal state Big Raman response of electron doped Fe 
SC [ini [21]. The Big quasi-elastic peak is a signature 
of the presence of critical fluctuations of the the nematic 
operator 

PBi, (q) = ^ (1) 

kcr 

whose susceptibility can be accessed in Raman scattering 
experiments in Big geometry [30]. Here ^Big (k) trans¬ 
forms as k‘^ — ky under the point group operations. Note 
that, Raman experiments thus imply that the nematic in¬ 
stability has a d-wave Pomeranchuk component. We take 
this as our starting point, and write a phenomenologi¬ 
cal interaction Hi = (q) (-q), with 

^ > 0, that can be tuned to study the Pomeranchuk- 
nematic QCP. The microscopic origin of g, which can 
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be due to spin or bond-current [m EiHsii or orbital 
[35l [ 37 H 39 ] or charge m fluctuations, is not relevant for 
the following discussion. Note that, fluctuations associ¬ 
ated with a Pomeranchuk instability are, by definition, 
intra-band terms. Furthermore, it is thought that the 
Big Raman response of electron-doped Fe SC is governed 
by the electron pockets mnn]. Consequently, to sim¬ 
plify the discussion, by we imply the fermions 

describing the electron pockets. With this understand¬ 
ing, we decouple l-Li by introducing the nematic bosonic 
Hubbard-Stratonovich field 0(q), and we get the action 

s = tJ 2 (xLm ( q , - 9) IV • • • • ( 2 ) 

The bare nematic susceptibility Xnemi^^'^^n) is the 
Fourier transform of (q, r) (q, 0))o, where 

(• • •)o implies average in the BCS ground state with 

p = 0. 

In Raman spectroscopy of Fe SC the momentum trans¬ 
fer by the photons is typically small compared to the fre¬ 
quency transfer {uj ^ where Ui? is typical Fermi ve¬ 
locity). Consequently, the Big Raman response function 
probes the imaginary part of the full nematic susceptibil- 
hy Xnem dynamical limit, and, in random 

phase approximation, is given by 

where we suppressed the momentum q = 0. 

For simplicity we evaluate the above at T = 0 and 
by ignoring lifetime broadening of the Bogoliubov quasi¬ 
particles. Later we comment about flnite-T and finite 
lifetime effects. Furthermore, we assume a constant su¬ 
perconducting gap A on the electron pockets (for other 
scenarii see [43l[44]). It is then straightforward to deter¬ 
mine the imaginary part of the bare Big Raman response 

m 


Imynem (^) 


^ LO^LU^ - i2Af 


( 4 ) 


This result is, up to the d—wave weighted electron-pocket 
density of states p = (k) (5 (ep - ek,eO> the 

usual one for a fully gapped superconductor [45]. The 
square root divergence for uj 2A from above is the 
well established pair-breaking peak in the electronic Ra¬ 
man spectrum of a superconductor [30]. For our consid¬ 
erations we also need the real part of Xnem (^) shown in 
Fig. ga). The important aspect of the real part occurs 
for |ci;| < 2A. In this regime it follows after Kramers- 
Kronig transformation 


^^Xnem 


(^) = X 


0 

nem,oo 


+ p 


(2A)^arcsin(^) 
uj^J (2A)^ — 


( 5 ) 



FIG. 1: (a) Real and imaginary parts of the bare BCS Raman 
response Xnem foi* an isotropic gap. The interacting response 
(Eq. HI in the main text) develops a resonance at Qr < 2A 
when the real part reaches the threshold 1/g (see Eq. [^ in 
text), (b) Development of the nematic resonance in the su¬ 
perconducting Big Raman response Rsig for different values 
of the dimensionless coupling constant A (see text). Note that 
for A = 0.3, 0.5, weak pair breaking continuum are visible. 


Besides the low-energy contribution, there is a weakly 
frequency dependent contribution from the high energy 
fermions, which we approximate by a constant Xnemoo- 
It can be absorbed in the dimensionless coupling A = 
pg/ (l — gXnem,oo)^ hcnccforth we take A as the tun¬ 
ing parameter that increases with increasing g. 

Because of its singular behavior, the leading contribu¬ 
tion to ReXnem {^) frequencies below 2A is the second 
term above. If we insert this result into Eqj^ we find a 
sharp pole in the renormalized Raman response once the 
condition 1 = ^fRex^gj^^ (fir) is fulfilled for fir < 2A, 
where the imaginary part vanishes. As illustrated in 
FigQa), the square root divergence of ReXnem (‘^) ^ 

approaches 2A from below, guarantees that the above 
condition can be satisfied and that a new exciton pole 
emerges. The resonance frequency of the nematic re¬ 
sponse is 


= 2Asin6>(A), (6) 

where 0{X) obeys the equation sin 2^ = 2X0. 

In the clean limit and at T = 0 iX, the mode is arbi¬ 
trarily sharp as it occurs below the particle-hole contin¬ 
uum that starts at 2A. Near the resonance the electronic 
Raman response is: 

(uj) = Zr6 {uj - Qr) (7) 

with spectral weight = ( (^) determined by ((x) = 
— /o ^ Zr must be compared with the 

ccvl——1) arcsin(cc) 

total weight 7r^pA/2 of the usual pair-breaking peak of 
the BCS theory. Near 2A the spectral weight vanishes 
linearly Since pA 1, the resonance 

mode soon acquires a weight comparable to the non¬ 
interacting BCS pair-breaking peak. A typical value for 
Qr = f 2A is 1.45. 

In Fig. Bb) we show the Raman response for different 
coupling constants. We included a small but finite in- 
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FIG. 2: Evolutions of the nematic resonance energy fir (a) 
and spectral weight Zr (b) with the dimensionless coupling 
constant A (dynamic limit). Because of hnite nemato-elastic 
coupling, Zr does not diverge, and Qr does not soften com¬ 
pletely, at Xqcp (see text). 


trinsic width F that may result from impurity scattering 
or thermal excitations of quasiparticles [44] . We see that 
for small nematic coupling constants the pair-breaking 
peak keeps its line shape but increases in weight. Once 
2A — fir is larger than F, a peak that is well separated 
from the continuum above 2 A emerges and sharpens with 
a Lorentzian line shape for larger coupling strength. The 
sharpness of the Raman resonance relies heavily on the 
opening of a full gap in the superconducting spectrum. 
Nodes of the gap around the electron pockets essentially 
wash out the nematic resonance [44] , making the presence 
of the nematic resonance a clear indication of the absence 
of nodes around the electron pockets. We also note that, 
already in the absence of any nematic correlations, dis¬ 
order induces a much weaker temperature dependence of 
the pair-breaking peak energy compared to the simple 
BCS result for A(T) [46], and we expect a similar effect 
for the T-dependence of the nematic resonance energy 

Lattice cutoff .— Next we discuss the fate of the ne¬ 
matic resonance as the system approaches the nematic 
QCP by increasing A. From Eq. Q it is tempting to de¬ 
duce that the QCP, defined by is given by 

A = 1. This, in turn, would imply that, using Eq. [^ the 
resonance frequency ^ 0, while the spectral weight 
of the resonance Zr ~ diverges. However, in prac¬ 
tice, this is not the case since, in addition to the purely 
electronic contribution of Eq. (i), there is a symmetry- 
allowed coupling to the lattice degrees of freedom given 
by m 


Hc = 't J df^rcj) (r) {dxu"^ - dyU^). (8) 

7 is a nemato-elastic coupling constant and u is the 
usual phonon displacement field. In the harmonic ap¬ 
proximation of the lattice, the effect of the coupling to 
the phonons is a renormalization of the nematic coupling 


constant 


+ (9) 

is the bare value of the orthorhombic elastic constant 

and q is a wave vector along the critical directions of the 

Brillouin zone [48]. Note that, the static limit ^fgtat = 
2 

^ (q, 0 ) = ^ ^ and the dynamic limit (q = 0, cj) = 

g do not commute [49]. Thus, the enhancement of the 
static nematic coupling constant due to the coupling to 
elastic degrees of freedom does not enter the dynamic 
Raman response which filters out the purely electronic 
contribution of the nematic response. 

On the other hand, it is the static limit that gov¬ 
erns thermodynamics and the actual nematic phase tran¬ 
sition. Thus, the condition for the nematic QCP is 
-^stat f 5 where Astat = P9sta.t/ (l - S'statXnem.oo) is the 
dimensionless nematic coupling constant in the static 
regime, renormalized by high energy excitations. Since 
^stat > if follows A < Agtat < 1- In other words, even 
at the nematic QCP fir does not soften to zero frequency, 
and concomitantly, Zr stays finite (see Fig. . In the 

symmetry-broken nematic phase (Agtat > 1) the behavior 
of {flr^Zr) is non-universal, and is beyond the scope of 
this work. 

Raman Experiments .— We now make contact with the 
experiments on Fe SC and in particular with electron 
doped BaFe 2 As 2 (Co-Bal22) [^ [191 EO]- We assume 
that decreasing Co doping x is equivalent to increas¬ 
ing A. In the superconducting state of Co-Ball2 the 
most salient feature of the Raman spectra is a peak ob¬ 
served in Big symmetry only [19]. Until now the Big 
peak has been attributed to a 2A BCS pair breaking 
peak coming from the electron pockets [DlpniiiT]. How¬ 
ever as first reported in Ref. [20], its Co doping depen- 
dance displays a striking behavior which is shown in Fig. 
I^a). Coming from the overdoped tetragonal side (x=0.1, 
Tc=20 K) of the phase diagram, the spectral weight of 
the peak increases dramatically upon approaching the or¬ 
thorhombic/nematic phase which occurs below x=0.065. 
It reaches a maximum at x=0.065 before collapsing in 
the orthorhombic phase at lower x. As stated above, the 
spectral weight of a simple BCS Raman pair-breaking 
peak is expected to scale as the gap energy A. How¬ 
ever the strong enhancement of around optimal dop¬ 
ing, where Tc is maximum, cannot be attributed to sub¬ 
stantial changes in A, since this would imply equivalent 
changes in Tc, whereas, in practice the Tc changes by less 
than 20% within the doping range considered. Moreover 
the peak energy in units of k^Tc actually softens upon 
approaching x=0.065, going from ^ 5 to ^ 4 /c^Tc (see 
Fig. [^c)), in a manner consistent with the nematic res¬ 
onance scenario. 

The above observations can be naturally explained by 
the nematic resonance scenario whose spectral weight Zr 
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FIG. 3: (a) Co doping (x) dependance of the Big Raman response in the superconducting state (SC) of Co-Bal22 [20]. The 
spectra have been subtracted by those just above Tc to highlight the SC induced features. x=0.065 corresponds to optimal 
doping (Tc=24.7 K) and the spectra in orange correspond to x compositions in the orthorhombic phase, (b) Phase diagram 
(T, x) of Co-Bal22 showing the evolution of Ts (orange) and of the electronic nematic Curie-Weiss temperature |To| (green) 
deduced from Raman experiments in the normal state m- Decreasing x content is equivalent as increasing the coupling A of 
our theory. Note that To is positive for x < Xc (filled symbols) and negative for x > Xc (open symbols), (c) Evolution of 
and Qr as a function of Co doping x in Co-Bal22. 


is strongly enhanced by the proximity of an nematic in¬ 
stability. In Co-Bal22 the position of the purely elec¬ 
tronic nematic instability Xc can be determined by look¬ 
ing at the enhancement of Raman Big nematic fluctua¬ 
tions observed in the normal state of the tetragonal phase 
m- The enhancement of the associated nematic suscep¬ 
tibility was found to follow a Curie-Weiss temperature 
dependence over a wide range of Co doping from x=0 
to x=0.1. As shown in Fig. [^b), the extracted Curie- 
Weiss temperature Tq is lower than Ts , the tetragonal to 
orthorhombic transition temperature, and extrapolates 
to T=0 K at Xc ^ 0.055. In a purely electronic model 
the spectral weight would diverge at Xc- However, 
as already stressed above, due to finite nemato-elastic 
coupling, the actual quantum critical point is moved to 
a higher doping xqcp, between x=0.06 {Ts=4:6 K) and 
x=0.065 {Ts=0 K), where and Clr are still finite as 
observed experimentally (Fig. |^b)) . 

We note that a similar divergence of the Big peak 
spectral weight has been observed recently in Co doped 
NaFeAs close to the boundary between the tetragonal 
and orthorhombic phase m- There the Big resonance 
was found to be almost Lorentzian below Tc, with a 
linewidth significantly smaller than in the case of Co- 
Bal22: 1 meV and 5 meV respectively near optimal 
doping. This difference possibly reflects a smaller cou¬ 
pling constant A in Co-Bal22 which would make the ne¬ 
matic resonance closer to the 2A continuum and there¬ 


fore broader. Recent ARPES data in Co-Bal22 indeed 
indicate that the Big peak energy is only slightly be¬ 
low 2A of one of the electron pocket [52]. Disorder may 
also play a role as optimally doped Co-Bal22 has a sig¬ 
nificantly higher Co content than Co-Nalll (see supple¬ 
mental materials for theoretical fits of the experimental 
lineshapes ISB). However, we stress that, while the line- 
shape of the resonance can be material dependent, the 
key feature of our theory, namely the enhancement of 
the spectral weight near the nematic QCP, is observed 
experimentally in both systems. 

The nematic resonance scenario captures remarkably 
well the salient features of the Big peak observed in the 
SC state of electron doped Fe SC: its nematic symmetry 
(Big) and the strong enhancement of its spectral weight 
near the nematic QCP. It also makes a direct link between 
the nematic response in the normal and superconducting 
state via Eq. close to Xc the enhanced nematic suscep¬ 
tibility of the normal state converts into a sharp nematic 
resonance at finite energy in the SC state because of the 
gap opening. 

Furthermore, our theory leads to the following three 
predictions which can be verified in future experiments, 
(i) From Eq. and assuming A to be T-independent, we 
expect that oc A upon varying temperature, (ii) Due 
to the nemato-elastic coupling of Eq. [^the orthorhombic 
elastic constant is renormalized to Cs = Cg = 

0). From Eq.j^and from Kramers-Kronig relation we de- 
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duce that the elastic softening 5Cs = Cg — scales as 
5Cs cx: upon approaching the nematic transition 

deep in the SC state, as a function of either doping or 
temperature, (iii) The nematic resonance will lead to 
single particle renormalizations (like peak-dip-hump fea¬ 
tures) everywhere on the electron pockets which should 
be visible in ARPES experiments. 

Conclusion .— In summary we have demonstrated the 
generic presence of a nematic resonance in the SC state 
of systems near a nematic quantum critical point like 
Fe SC. The very existence of this resonance relies only 
on the presence of a fully gapped Fermi pocket and the 
proximity of a nematic quantum critical point. Existing 
Raman data in the Big symmetry channel of several Fe 
SC systems like Co-Bal22 and Co-Nalll indicate that 
this resonance actually dominates the Raman response in 
the SC state. It provides the most striking manifestation 
of nematicity in the superconducting properties of Fe SC. 
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